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Abstract 

We define a multivariate polynomial that generalizes several interlace 
polynomials defined by Arratia, Bollobas and Sorkin on the one hand, 
and Aigner and van der Hoist on the other. We follow the route traced 
by Sokal, who defined a multivariate generalization of Tutte's polynomial. 
We also show that bounded portions of our interlace polynomial can be 
evaluated in polynomial time for graphs of bounded clique-width. Its full 
evaluation is necessarily exponential just because of the size of the result. 

Keywords : Interlace polynomial, multivariate polynomial, monadic second- 
order logic, clique- width. 
Support^ 



1 Introduction 



Many polynomials associated with graphs, matroids or combinatorial maps 
count configurations in these objects. We take here "configuration" in a wide 
sense. Typical examples are colorings, matchings, stable subsets, subgraphs. In 
many cases, a weight is associated with the considered configurations : number 
of colors, cardinality, number of connected components or rank of an associated 
subgraph. 

A multivariate polynomial, as the one defined by A. Sokal, that generalizes 
Tutte's two variable polynomial, not only counts configurations, but also enu- 
merates them, together with their associated weights. Multivariate polynomials 
may have recursive definitions, the specializations of which give the classical 
recursive definitions for the corresponding polynomials. (Specialization is the 
inverse of generalization, these notions are explained below) . 

We think that a recursive definition at the multivariate level makes better 
understand what is going on at the usual level of one or two variable polynomials. 



^This work has been supported by the ANR project GRAAL and by a temporary position 
of CNRS researcher. 
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Furthermore, multivariate polynomials may have static definitions based on 
properties of configurations and weights expressible in second-order logic or 
better in monadic second- order logic. Why better ? Because this yields fixed 
parameter tractable algorithms, where the parameters arc the tree-width or the 
clique-width of a graph, or the branch-width of a matroid. This consequence 
has been explained in [CMR, Mak04, Mak05]. We apply these ideas to the 
interlace polynomials defined first under this name in [ABS], then generalized 
in [ABS04b] and [AvH] (and previously defined under another name, see Las 
Vergnas [LV]), that is, we give a common multivariate generalization and show 
that is is definable by monadic second-order formulas. 

What do we mean by "generalize" ? If two polynomials P{G) and Q{G) are 
associated with each graph G of a certain type, we say that P{G) is more general 
that Q{G) and that Q{G) is a specialization of P{G) if Q{G) can be obtained 
from P{G) by a substitution of fixed polynomials to the variables of P{G); 
such a substitution may be also combined with multiplication by polynomials 
depending, say, on the number of vertices and/or edges. We do not try to 
propose here a most general definition. Applications of the idea will suffice. 

A multivariate polynomial is one with indeterminates depending on the ver- 
tices or the edges of the considered graph (such indeterminates are sometimes 
called "weights", because they make possible to evaluate the polynomial with 
distinct values associated with distinct vertices or edges). Sokal's multivariate 
Tutte polynomial of a graph G = (V, E) is defined by : 

Z{G)^EACE^''^''^^^^Ue^AVe 

where G[A] is the subgraph of G with set of vertices V and set of edges A, 
k{G[A]) is the number of its connected components. This polynomial belongs 
to Z[u,Ve ; e e E]. An indeterminate Ve is associated with each edge e. The 
indeterminates commute, the order of enumeration over each set A is irrelevant. 
For two graphs G and G' with sets of edges in bijection, we have Z{G) = Z{G') 
(where the variables indexed by edges of G and G' correspond via the considered 
bijection) iff | V(G) |=| V{G') \ and their cycle matroids are isomorphic (via 
the same bijection between edges). This observation explains what information 
is contained in this polynomial about the considered graph. 

The polynomial Z{G) is more general than Tutte's two variable polynomial 
T(G, X, y) because (see [Sok] for details) : 

T(G, X, y) = {{x - lf^^\y - l)\^\)-^a{Z{G)) 

where a is the siibstitution : 

[u := {x — l){y — l);Ve := y — 1 for all e e E]. 

Conversely, one can express Z'(G), defined as Z(G) where every indetermi- 
nate Ve replaced by the same indeterminate v, in terms of T(G, x, y) in a similar 
way. Hence, Z'(G) and T(G) are equivalent both in expressive power and for 
the complexity of their computations. 
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In this article, wc define a multivariate polynomial, that generalizes the 
different interlace polynomials defined for graphs in [ABS] , [ABS04b] and [AvH] , 
and also, the independence polynomial surveyed in [LM]. In our polynomial, 
a configuration is a set A of vertices and its associated value is the rank of 
the induced subgraph. By rank we mean, as in these articles, the rank of 
the adjacency matrix with respect to the field GF(2). Actually, in order to 
generalize a polynomial defined in [AvH], we introduce a second argument in a 
configuration, namely a set B of vertices corresponding to "toggled loops" : if 
a vertex in B has no loop we add one ; if it has a loop, we delete it. Then we 
evaluate ranks. 

We find for this polynomial a recursive definition, somewhat more compli- 
cated than the usual ones based on contracting and deleting edges. The poly- 
nomials of [ABS], [ABS04b] and [AvH] are specializations of ours, and we find 
their recursive definitions as the corresponding specializations, sometimes with 
the necessity of proving nontrivial properties. 

This approach is based on static definitions from which recursive definitions 
can be constructed, and not the other way around. Here we follow Sokal who 
considers recursive definitions as secimdary. It is dcvelopped in a fundamental 
perspective in [CGM]. Let us say to have a representative picture of the most 
general case that the polynomials we obtain in this way are of the form : 

where C ranges over all configurations of a multiset T{G), nc is the number 
of occurrences of C in r(G), vc is a monomial (like OeGA above 
polynomial Z{G)) that describes configuration C, and f{C) is the weight of 
C. Polynomials of this form have necessarily positive coeffficients. We are 
especially interested in cases where T{G) and / can be expressed by monadic 
second-order formulas because of algorithmic applications developped in [OMR, 
Mak04, Mak05]. 

Such polynomials are usually of exponential size. We may ask for a polyno- 
mial evaluation of the truncation of P{G) i.e., of its part limited to the mono- 
mials associated with configurations of a given size. For such purpose, having 
a description of T{G) and of the weight function / by monadic second order 
formulas is helpful, as we will see in the last section, because this yields fixed pa- 
rameter tractable algorithms, where the parameter is tree-width or clique-width. 

Summary of the article : 2. Definitions and basic facts, 3. A multivariate 
interlace polynomial and its recursive definition, 4. Specializations to known 
polynomials, 5. Polynomial time evaluation for graphs of bounded clique- width, 
6. Conclusion, 7. References, 8. Appendix : The Tutte polynomial 
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2 Definitions and basic facts 



Graphs are finite, simple, undirected, possibly with loops. A graph is defined as 
a pair G — {Vq, Aq) of a set of vertices Vq and a symmetric adjacency matrix 
Aq over GF(2). We omit the subscripts whenever possible without ambiguity. 
The rank rk{G) of G = {V, A) is defined as the rank rk{A) of A over GF(2); its 
corank (or nullity) is n(G) := n{A) :=[ V \ —rk{A). The empty graph has 
rank and corank 0. 

The set of looped vertices of G (the vertices i such that A{i, i) — 1) is denoted 
by Loops{G). For a in V, we let N{G,a) be the set of neighours b of a, with 
h ^ a. (A looped vertex is not a neighbour of itself). 

If X is a set of vertices, we let G — X denote G[V — X], the induced subgraph 
of G with set of vertices V — X. 

We denote by GVX the graph obtained by "toggling" the loops in X, i.e., 
Vgvx ■■= Vg and : 

AGvx{iJ):=l-AG{iJ) if i = j G X, 
AGVx{i,j)-=AG{i,j) otherwise. 

We write G = H (B K ii G is the union of disjoint subgraphs H and K. 

For two graphs G and H we write H = h{G) and we say that they are 
isomorphic by his his a. bijection of Vq onto Vh and AH{h{i), h{j)) = ^^(i, j) 
for all i and j. 

Pivoting and local complementation 



We recall the precise definitions of operations like local complementation 
and pivoting, because there are some variants in articles. 

For a and b distinct vertices of G we define the graph H = G""^ as follows : 
Vh -—Vg and 

Aniijj) := 1 — AG{i,j) if the following holds : 
{hj} ^ {^1 b} = and 

{either i G A^(G, a) - N{G, b) and j G iV(G, b), 

or j G N{G,a) - N{G,b) and i G N{G,b), 

or iGN{G,b)~N{G, a) and jGNiG,a), 

or j e iV(G, b) - N{G, a) and i G N{G, a)}. 
In all other cases, we let AH{i,j) ■= AG{i,j)- 

This operation is called pivoting on a, b. It does not depend on whether a 
and b are loops or are adjacent. 

Local complementation : for a vertex a of G we define = G" as follows : 

Vh ■= Vg and : 

AH{i,j) := 1 — AG{i,j) if i,j G N{G, a), including the case i = j. 
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AH{i,j)-=AG{i,j) otherwise. 
Another notion of local complementation is defined by : 

G*a= {GVN{G, a))" G"ViV(G, a). 

It "toggles" the (non loop) edges of G[N{G, a)]. It is used for graphs without 
loops in the characterization of circle graphs and in the definition of vertex- 
minors. ([Bou], [Oum], [CouOum]). 

We write a — 6 to express that a and b are adjacent both without loops, and 
— 6 to express the same with a looped and 6 not looped, and — b^^ if a and 
b are looped. These operations satisfy properties listed in the following lemma : 

Lemma 1 : For every graph G = {V,A), for distinct vertices a,b and all 

sets of vertices X, Y wc have : 

(1) (G")° = G; G"** = G''" ; (G"'')"'' = G ; 

(2) G"'' /i(((G")'')"Va) if 6 e N{G, a) and h is the permutation of V that 
exchanges a and 6 ; G°'^ - a - b = ((G")'')'^ - a - b. 

(3) (G"*')*' = /i((G")''Va) if a, 6, are as in (2) ; {G"'')'' -a-b = {G")'' -a-b. 

(4) GVX'"'' = G°^VX ; GVX" = G^VX ; GVX[Y] = G[Y]V{X nV). 

(5) G[X]»'' = G''^[X] ■ G[X]'' = G"[X] if a and b are not in X. 

Proof : (1), (4), (5) are clear from the definitions. 

(2) is a well-known fact about pivoting and local complementation. See for 
instance [ABS04b]. 

(3) This is a consequence of (1), (2) and (4) : 
(G"'')^ = (/i(((G»)'')"Va))*' 

= /i(((G«)^)«Va)'') = /i((((G«)'')«)«Va) = /i((G")''Va).n 

Computating ranks of graphs. 

Lemma 2 : For every graph G, for distinct vertices a, 6 we have : 

(1) rfc(G) = 1 + rA:(G" - a) if a G Loops{G) ; 

(2) rfc(G) = 2 + rA;(G"^ - o - 6) if a - 6 ; 

(3) rfc(G - a) = rfc(G"'' - a) if a - &; 

(4) rfc(G) = 2 + rA:((G")^ ~ a - b) = 1 + rk{G''^ -b)iia^-b . 
Proof : (l)-(3) arc proved in [ABS04b]. 

(4) We note that {G^f - a - b = {G" - a)'' - b and that {G" - a) has a 
loop on b. Hence by using (1) twice : 

rkiiG^f -a-b)= rk{{G'' - a)^ - b) = rfc(G" - a) - 1 = rk{G) - 2. 
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For the second equality : 

rk{G) = rk{{G^Y - a - 6) + 2 = rk{{{{G^fff - a - 6) + 2 

= rk{{{{G^Yf -af-b) + 2 = rk{{{G''Yf - a) + 1 

because {{G^YY ~ ^ ^^'^ ^ loo^i on 6, hence : 

rk{G) = rk{h{G^''Vb) - a) + 1 = rk{h{G''''Vb -b)) + l 

= rk{G'"'Vb -b) + l= rk{G"-^ - 6) + 1 

because G^'^Vb -b = G°-^ - b.U 



3 The multivariate interlace polynomial 



Polynomials have integer coefficients. Following Sokal [Sok] wc call multivariate 
polynomials those with indeterminates Xa, Da^ Za,--- associated with vertices a 
of the considered graph G. We will denote by Xg the set of such indeterminates 
for X = {x,y,z, ...}. They are the G-indexed indeterminates. We denote by U 
a set {u,v,w, ...} of "ordinary" indeterminates not associated with elements of 
graphs. 

By a polynomial P{G), we mean a mapping P that associates with a graph 
G a polynomial in Z[[/ U Xg] such that if h is an isomorphism of G onto H, 
then P{H) is obtained from P{G) by the substitution that replaces Xa by a;^(a) 
for every Xa in X^. 

A specializing substitution is a substitution that replaces an indeterminate 
from a finite set U = {u,v,iu, ...} by a polynomial in Z[{7], and a G-indexed 
indeterminate Xa in Xg, by a polynomial in Z[t/ U {ya \ y G X}], the same for 
each a. For an example, such a substitution can replace Xa by ya{x — 1)^ — 
3ZaU + 1, for every vertex a of every graph. If cr is a specializing substitution, 
then a o P, defined by a o P{G) = a{P{G)) is in this sense a polynomial. 

For a set A of vertices we let xa abbreviate the product (in any order) of 

the commutative indeterminates Xa, for a in A. If A = 0, then = 1. If i? is 
a set of subsets of G, then the polynomial J2AeB -^^ describes exactly B. If B 
is a multiset of sets, then it is described by J^AeB '^l^) ' where n{A) is the 
number of occurrences oi A in B. 

Definition 3 : The multivariate interlace polynomial. 

For a graph G we define 



where A, B range over subsets of V. Hence B{G) € Z[{u, v} U Xg ] where 
X = {x,y}. 
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The interlace polynomial q of of [ABS04b] is obtained from B{G) by a sub- 
stitution: q(G;x,y) = a{B(G)) where a is the substitution : 

[u := X — l;v := y — l;Xa ■= l,ya := for all a gV], 

and the polynomial Q of [AvH] , defined for graphs without loops is Q{G, x) = 
t{B{G)) where r is the substitution 

[u := l;v := X — 2; Xa ■= Va '■= 1 for all a gV]. 

These polynomials are actually defined recursively in these articles (which 
raises a problem of well-definedness), and then proved to be equal to the poly- 
nomials (t{B{G)) and t{B{G)). 

For more clarity with variable names, we will write q{G) and Q{G) with 
variables u' and v' instead of x and y. Hence q{G;u',v') = a{B{G)) where a 
is the substitution : 

[u := u' — l;v := v' — 1; Xa ■= 1, ya '■= for all a e V], 

and Q{G,v') = t{B{G)) where r is the substitution 

[u := l;v := v' — 2;Xa '■= ya '■= 1 for all a G V]. 

Let Bi{G) be the polynomial obtained from B{G) by replacing w by 1. 

Lemma 4 : For every graph G and every set T of vertices : 

(1) B{G) = eiBiiG)) where : 

6 := [u := u'ir^;xa := vxa, Va vya for all a G V], 

(2) B{GVT) = i^{B{G)) where : 

/X := [Xa ■■= ya, ya ■= Xa for all a G T]. 

Proof : (1) Clear. 

(2) We observe that GWTWB[A U B] = GW{A' U B')[A U B] where A' = 
AnT,B' = B - BnT. The result foUows.D 

We will write : B = 6 o Bi. The polynomial B{G) can thus be "recovered" 
from Bi{G). Since every graph G is GiVT for some T with Gi without loops, 
we have B{G) = ijl{B{Gi)) where /U is as in Lemma 4. Hence, it is enough to 
know B{G) for graphs G without loops. However, the recursive definitions to 
be considered below will introduce graphs with loops in the recursive calls. 

Properties of polynomials 

The polynomials q and Q defined above satisfy some properties for all graphs 

G : 

q{G -a)-q{G-a-b) = q{G''^ - a) - ^(G"'' - a - 6) if a- 6(1) 
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Q{G*a)=Q{G) (2) 
Q{&^^) = Q{G). (3) 
Do these equalities hold for B{G) ? The answer is no for (2) and (3) as a 

consequence of the next proposition (and also for (1), see below Counter-example 

14). 

Proposition 5 : A graph G and its polynomial B{G) can be reconstructed 

from p{B{G)) where p := [v := l;ya ■= for all a G V]. 

Proof : For every set of vertices A, the rank of G[A] is the unique integer 
n such that Xyiw" is a monomial of p{B{G)). Now a vertex a has a loop if 
rk{G[a]) ~ 1, and no loop if rk{G[a\) = 0. Hence, we obtain Loops{G) from 
p{B{G)). Using this information, we can reconstruct edges. 

If a and b are not looped, they are adjacent iff rfc(G'[{a, h}] = 2, otherwise 
rfc(G[{a,6}] =0. 

If one of a, 6 is looped, they are adjacent iff rk(G[{a,b}] = 2, otherwise 
rfc(G[{a,6}] = 1. 

If both are looped, they are adjacent iff rk{G[{a, b}] = 1, otherwise rk{G[{a, 
b}] = 2.n 

This proof shows how strange is the behaviour of the rank function on graphs. 
It follows that identities (2) and (3) cannot hold for B and even for po B. 

Question : By which algebraic transformations can B{G) be expressed in 
terms of p{B{G)), in a uniform way, valid for all graphs G ? 

3.1 Recursive definition 



We now determine a recursive definition, (also called a set of reduction formulas) 
of B(G) from which will follow the recursive definitions of [ABS] and [AvH]. 

We let a denote the graph with one non-looped vertex a, and the similar 
graph with looped vertex a. 

Lemma 6 : For every graph G, for every graph H disjoint from G we have: 

(1) Bi0) = 1 

(2) B{G ®H) = B{G) ■ B{H) 

(3) B{a) = 1 + XaV + DaU 

(4) B(a^) = l + XaU + VaV. 

Proof : Easy verification from the definitions. □ 

The more complicated task consists now in expressing B{G) in the case where 
a and b are adjacent (this is necessary if no rule of Lemma 6 is applicable). We 
will distinguish three cases : a — b, — b, and — b^. 
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For a graph G and disjoint sets of vertices A and B, wc let m(G, B) denote 
the monomial xaVbu'^^^^^'^^'^^^'iv^'^^^'^'^'^^^^ so that B{G) is nothing but the 
sum of them over all pairs A, B (the condition An B = will be assumed for 
each use of the notation m(G, A, B)). 

For distinct vertices a, b, two disjoint sets A, B can contain a, b or not ac- 
cording to 9 cases. We let i e {0, 1, 2} mean that a vertex is in F — (A U B), in 
A or in i? respectively. Let Bij be the sum of monomials m{G^ A, B) such that 
i tells where is a, and j tells where is 6. For an example : Bq2 is the sum of 
monomials m{G, A, B) such that a & V - {A\J B), b & B. 

Claim 7 : Let G with a — b. 

(1) Boo = B{G-a-b) 

(2) Bii = XaXbu'^ ■ B{G<'^ -a-b). 

(3) B20 = Vau ■ BiG" -a-b); B02 = ytu ■ B{G'> -a-b); 

(4) S12 = XaVbU^ ■ B{{G^)'' -a-b); B21 = XbVau'' ■ -a-b). 

Proof : (1) Clear from the definitions. 
(2) A monomial of Bn is of the form : 

m(G, A, B) = a:^y5u-fc(GVB[^uB])^n(Gvs[AuB]) 
with a,b & A (because of the subscript 11). By Lemma 2(2) we have : 
rk{GS/B[A U S]) = 2 + rk{GVB[A U Bf^ -a-b). 

But GWB[A U BY'' -a-b= (G"*" -a- b)WB[A' U B] where A' = A- a-b 
(we use here Lemma 1(4,5)). Hence : 



because the set of pairs A',B C V — a — b such that A' and B are disjoint 
coincides with the set of pairs {A — a — b),B such that A,BCV,A and B are 
disjoint subsets of V and a,b G A. 

(3) The proof is similar. A monomial of B20 is of the form (1) above with 
a G B, b ^ Au B (because of the subscript 20). By Lemma 2(1) we have: 



m(G, A, B) = XaXbU^ ■ m{G"^ -a-b, A', 



B). 



It follows that 



Bn = XaXbU^ ■ BiG"'' -a-b) 



rk{GVB[A U B]) = 1 + rk{GVB[A U B^ - a) 



because a is looped in GVB[A U B]. But : 



GVB[AUB] 



a 



— a = 



{GS/a" -a-b)VB'[AuB'] 
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because b ^ AU B, where B' = B — a. (By Lemma 1). Clearly, (GVa" - 
a-b) = -a-b). Hence m(G, A, B) = y^u ■ m{G'' -a-b,A, B'). It follows 
that : 

B20 = Vau ■ B(G" -a-b) 

because the set of pairs A,B' C V — a — h such that A and B' are disjoint 
coincides with the set of pairs A, {B — a) such that A,BCV,A and B are 
disjoint subsets of V,a G B and b ^ Au B. The case of B02 is obtained by 
exchanging a and 6. 

(4) A monomial of B12 is of the form (1) above with a € A, b € B. By 
Lemma 2(4) we have : 

rk{GVB[A U B]) = 2 + rk{{GWB[A U B]^)" -a-b) 

because 6^ - a in GVB[A US]. We have : 

{GVB[A U B]'')" -a-b= ((G'')" -a-b) VB'[A' U B'] 

where A' = A - a, B' = B - b. Hence : 

m(G, A, B) = XaVbu'^ ■ m((G'')" -a-b, A', B'). 

It follows that : 

B12 = • B((G'')« - a - 6) 

because the set of pairs A' ,B' <ZV — a-b such that A' and B' are disjoint 
coincides with the set of pairs {A — a),{B — b) such that A, B C V , A and B are 
disjoint a £ A and b £ B. The case of -B21 is obtained similarily by exchanging 
a and b.O 

The next claim establishes linear relations between some polynomials Bij . 

Claim 8 : Let G with a — b. 

(1) B{G - a) = Boo + Boi + B02 

(2) B(G - 6) - Boo + Bio + B20 

(3) uya ■ B{G'' - a) - B20 + B21 + B22 

(4) uyi, ■ B{G'' -b)= B02 + B12 + B22 

Proof : (1), (2) Clear from the definitions. 

(3) From the definitions, B20+B21+B22 is the sum of monomials rn(G, A, B) 
such that a G B. We have : 

rA;(GVB[^ U B]) = 1 + rk{GVB[A U B]" - a) 



10 



by Lemma 2(1). But : 

GVB[A yjBf -a = {{GVaf - a)VB'[A U B'] (where B' = B - a) 
= {G"- -a)VB'[AuB']. 

This gives the result with the usual argument. 
(4) Similar to (3) by exchanging a and b.D 

If we collect the equalities of Claims 7 and 8 we have 10 definitions or linear 
equalities for 9 "unknowns". This is enough for obtaining B{G). We get thus : 

B{G) = (i?oo + -Bio + -B20) + {-Boi + Bii + B21} + (i?02 + -B12 + B22) 
= B{G -b) + {Bqi + XaXbu' ■ BiG"'' -a-b)+ XbVaU^ ■ B{{G''f -a- b)}+ 
+ybU-B{G^ -b). 

Then Boi = B{G - a) - Bqq - B02 

= B{G -a)- B{G -a-b)-ybu- B{G^ -a-b) . 

We obtain, after reorganization of the expression and a bit of factorization ; 

Lemma 9 : Let G with a — b. Wc have : 
S(G) = XbU^Xa ■ B{G-b -a-b)+ya- S((G")'' -a- b)}+ 
+ybu{B{G'' - 6) - B{G'' -a- b)}+ 
+B{G -a)+ B{G - 6) - B{G -a-b). 

Considering B22 for which we have two expressions, we get : 

Corollary 10 : Let G with a — b. 

yb{B{G'> -b)- B{G^ -a-b)-XaU- B{{G°'f -a-b)} 

= yaiBiG" -a)- B{G"- -a-b))- XbU ■ B{{G^)'' -a- b)}. 

Next we consider the cases a^ — b and — b^. Actually, Lemma 4(2) will 
shorten the computations. 

Lemma 11 : (1) Let G with a — b^. 

B{G) = ybu^xa ■ B{G-b -a-b) + Va- B{{G'')^ -a-b)} 
+Xbu{B{G'' -b)- B{G'' -a-b)} 
+B{G -a)+ B{G -b)- B{G -a-b). 

(2) Let G with - 5^. 

B{G) = ybu'iya ■ B(G'^'' -a-b) + Xa- BUG'^)" -a-b)} 
+Xbu{B{G'' - 6) - B{G'' -a-b)} 
+B{G -a)+ B{G -b)- B{G -a-b). 

Proof : (1) We have G = G1V6, Gi = GV6, where in Gi we have a — 6 so 

that Lemma 9 is applicable. 

We get then, letting (3 be the substitution that exchanges Xb and yb : 

B(G) = /3(B(Gi)) = 



11 



= Vbu'^ixa ■ BiGVb"'' -a-b)+ya- BUCVb^f -a- b)}+ 
+Xbu{B{GVb'' - 6) - B{G\/b'' -a- b)}+ 
+(3{B{GVb - a)) + B{GVb -b)- B(GV6 - a - 6) 

= VbU^Xa ■ BiG'^" -a-b) + ya- B{{G'^f - a - b)] 
+Xbu{B{G^ -b)- B{G^ -a- b)} 
+B{G -a)+ B{G -b)- B{G -a-b), 

because we have GVb"^ - a-b = G"'' - a-b and B{GVb''^ -a-b) has no 
occurrence of an indeterminate indexed by 6, 

because we have (GSIb^f -a-b= {G^f - a - 6 and B{{GS7b''f -a-b) 
has no occurrence of an indeterminate indexed by 6, 

and similar remarks apply to GVb^ — b, to G\7b^ — a — 6, to GVb — b, and 
to GVb - a- b. Finally, we have (3{B{G'\/b — a)) = B{G — a) hy Lemma 1 and 
Lemma 4. 

(2) Very similar argument. □ 

We can now sum up the results of lemmas 6, 9, 11 into the following propo- 
sition, where the three cases arc collected into a single one with help of the little 
trick of introducing "metavariables" ZcWc for each c S F : 

Zc = Xc and Wc = Uc if c is not a loop, 

Zc = Vc and Wc = a;c if c is a loop. 

Proposition 12 : For every graph G, for every graph H disjoint from G, 
every vertex a, we have : 

(1) B{0) = 1 

(2) B{G (B H) = B{G) ■ B{H) 

(3) B{a) = 1 + XaV + UaU 

(4) B{a^) = l+XaU + VaV 

(5) B{G) = Zbu^iza ■ BiG"'' -a-b) + wa- B{{G"f -a-b)} 

+Wbu{B{G^ - 6) - B{G^ -a-b)] 
+B{G -a)+ B{G - 6) - B{G -a-b). 

if 6 G N{G,a). 

Proof : Immediate consequence of Lemmas 6,9,11. □ 
We have an even shorter definition : 

Corollary 13 : For every graph G every vertex a, we have : 

(1) B(0) = 1 

(2) B{G) = (1 + ZaV + Wau)B{G - a) if N{G, a) = 0, 

(3) B{G) = Zbu'^{zaB{G''^ -a-b)+ WaB{{G'')^ -a-b)} 

+Wbu{B{G'' -b)- B{G'' -a-b)} 
+B{G -a)+ B{G -b)- B{G -a-b). 

if 6e Ar(G, a). 
Counter-example 14 : 

It is proved in Proposition 8 of [ABS04b] that if a — 6 in G then : 
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q{G -a)-q{G-a-b) = - a) - -a-b) 



This is not true for B. Take G = c — a — b — d. Note that G"'' is G augmented 
with c — d. Assume we would have : 

B{G -a)~ B{G -a-b)^ B{G'''' - a) - B(G"'' -a-b). (*) 

In the left handside, we have a single monomial of the form UbUcXdu" for 
some n, and it must be from B{G — a) because b is not in G — a — 6. This 
monomial is ybVcXdU^ because rk{c^ ffi (d — b^)) — 3. In the right handside we 
have the monomial ybVcXdU^ because rk{c^ — d — b^) = 2. Hence we cannot have 
Equality (*).□ 

In such a case, we can ask what is the less specialized (or most general) 
substitution a such that the corresponding equality is true for a o B ? Some 
answers will be given below. We prove actually a more complicated identity. 

Proposition 15 : If a — 6 in G then : 

B{G - a) - B{G -a-b)- BiG"^ - a) + B{G'''' - a - b) = 
ybu{B{G^ -a-b)- B{{G''f -a- b)}. 

Proof : We use the notation and some facts from Claims 7 and 8 : 

B{G -a)- B{G - a - 6) = Boi + -B02 = -Boi + VbU ■ B{G'> -a-b). 

We let Bqi ^iid Bq2 denote the polynomials Bqi and B02 relative to (G"**, a, b) 
instead of to (G, a, b). Then we have : 

BCG"*" - a) - B(G«^ - a - 6) = B§^ + B§l = Bg^ + ybU ■ B{{G''^f -a-b). 

We have by Lemma 1 : {G"'')^' -a-b= {G'^f -a-b. 
On the other hand, Bqi the sum of monomials : 

m(G«^ = a;^y3„'-fe(G»^VS[AUB])^n(G»^VS[AUB]) 

for disjoint sets A, B such that a^A\JB,b<EA. But for such A, B : 

G"''VB[A UB] = GWBIA UBU a]"'' - a. 

Hence, using Lemma 1 and Lemma 2(3) : 

rfc(G"''VB[A U B]) = rk{GVB[A UBU a]"'' - a) 
= rk{GVB[AuBUa]-a) 
= rk{G\/B[AuB]). 

We have also n(G"''VS[A U B]) = n{GVB[A U B]). Hence, m(G"^ A, B) = 
m{G, A, B) and B^'l = Bqi. 

13 



Collecting these remarks we get : 

B{G - a) - B{G -a-b)- B{G''^ - a) + B{G°-^ ~a-h) 
= Bo2 - -Bol 

= VbU ■ B{G^ -a-b)-ybu- B{{G°'f -a- h)n 

We note for later use that Identity (*) of Counter-example 14 holds if either 
w = or 7/5 = for all h. 

A polynomial P in Z[X] is said to be positive if the coefficients of its mono- 
mials are positive. A polynomial in the sense of a mapping from graphs to 
polynomials is positive if P{G) is positive for every graph G. It is clear from 
Definition 3 that B is positive. This not immediate from the recursive defini- 
tion of Corollary 13 because of two substractions in the right handside of (3). 
However, one can derive from Corollary 13 a stronger statement that is not 
immediate from Definition 3. 

Proposition 15 a : For every graph G and every vertex a, the polynomials 
B{G) and B{G) - B{G - a) are positive. 

Proof : By induction on the number of vertices of G, one proves simulta- 
neously these two assertions by using Corollary 13. 
In case (2) we have: 

B{G) — B{G — o) = (1 -I- ZaV + Wau)B{G — a) and in case (3) we have 
B{G) - B{G -a)= Zbu^{za ■ B{G''^ -a-b)+Wa- B{{G''f -a-b)] 
+Wbu{B{G^ -b)- B{G^ -h-a)} 
+B{G -b) - B{G -b-a), 
which gives with the induction hypothesis that B{G) — B{G — a) is positive. 
So is B{G) since, again by induction, B{G — a) is positive. □ 



4 Specializations to known polynomials 



We have already observed that the polynomial q of [ABS04b] is q{G\ u', v') = 
a{B{G)) where a is the substitution : 

[u := u' — l;v := v' — 1; Xa ■= 1, Va '■= for all a €V], 

and that the polynomial Q of [AvH], defined for graphs without loops is 
Q{G,v') = t{B{G)) where r is the substitution 

[u := 1; V := v' — 2; Xa ■= Va '■= 1 for all a G V]. 

Both are actually specializations of the following two. We let : 

By^Q{G) := ao{B{G)) where ctq is the substitution := for all a £ V], 
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and 

Bx^y{G) := a^{B(G)) where a= is the substitution [ya := Xa for all a G F]. 

Polynomials B , Bx=y, By^Q arc by definition positive. This is also the case 
of Q) but this is not obvious from the above definitions. This fact follows for 
Q from the recursive definition established in [AvH] that we will reprove in a 
different way, but it does not from the definitions given in [ABS04b].(In the 
appendix we give a proof that a multivariate version of the Tutte polynomial is 
positive.) 



4.1 Fixed loops 

The polynomial By^o{G) corresponds to "fixed loops" : it does not describe 
what happens when some loops are "toggled". It can be rewritten : 

Cleanly q{G;u',v') = a' {By^Q{G)) where a' is the substitution : 

[u := u' — 1; V := v' — 1; .x„ := 1 for all a G V]. 

Proposition 16 : For every graph G every vertex a, we have : 

(1) i?,=o(0) = 1 

(2) By=o{G) = (1 + Xav)By=o{G — a) if N{G, a) = and a is not a loop, 

(3) By^aiG) = XaU ■ i?y=o(G° — a) + By=o{G — a) if a is a loop, isolated or 
not, 

(4) By=o(G) = xbxau^ ■ Sj,=o(G'^' -a-b)+ 

+By=o{G -a)+ By=o{G -h)- By=o{G -a-b) if a - 6. 

Proof : (1), (2), (4) : Immediate from Corollary 13. 

(3) If a is isolated, this follows from Corollary 13 (2). Otherwise, using 
the notation of the proof of Claim 7 we note that By=o{G) is the sum of 
monomials m{G, A, 0) ; those such that a ^ A yield B{G — a), the others yield 
XaU ■ By^o{G'^ — a) since : 

rk{G[A]) = rfc(G[A]» - a) + 1 = rfc((G» - a)[A - a]) + 1 

by Lemma 2(1). This gives the result, however, it is interesting to see what 
gives Lemma 11. The two cases — b and — b^ yield the same equality. 
By^oiG) = XaulBy^oiG" -a)- By^oiG'' -a-b)} 

+By=oiG -a)+ By=o{G -b)- By=o{G -a-b). 
Hence we have to check that : 

XaU ■ By^oiG" -a-b) = By=o{G - 6) - Sj,=o(G -a-b). 
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This is nothing but Assertion (3) applied to H = G — b. Hence (3) can be 
estabhshed by induction on the size of G, with help of Lemma 11, and without 
repeating the analysis of the monomials m{G, A, 0).D 

This proposition yields, with easy transformations, the following recursive 
definition of q : 

(ql) q{G) = w'" if G consists of n isolated non-looped vertices, 

(q2) q{G) = [u' — l)q{G°' — a) + q{G — a) if a is a loop, isolated or not, 

(qS) q{G) = [u' - lfq{G''*> -a-b)+ 

+q{G -a) + q{G-b)-q{G-a-b) if a - 6 

which is not the one of [ABS04b] Proposition 6. The definition of this article 
replaces (q3) by : 

(q3') q{G) = {{u' - 1)^ - 1)^(6'"'' -a-b)+ 

+q{G -a) + qiG"'' - b) if a - 6 

However, we have the following corollary of Proposition 15 that generalizes 
Proposition 8 of [ABS04b] : 

Corollary 17 : If a — 6 in G then : 

By=o{G - a) - By=o{G -a-b) = By^oiG'^^ - a) - By^oiG^^^ -a-b). 



Proof : Immediate from Proposition 15 since 2/6 = for all b. □ 
We get thus the following corollary, and (q3') is equivalent to (q3). 

Corollary 18 : For every graph G every vertex a, we have : 

(1) i?,=o(0) = 1 

(2) By=o{G) = (1 + Xav)By^o{G — a) if N{G, a) = and a is not a loop, 

(3) By^o{G) = XaU ■ By^Q{G"' —a) + By^o{G — a) if a is a loop, isolated or 

not, 

(4) By=o{G) = {x^xau" - l)Sj,=o(G"'' - a - &)+ 

+By=o (G - o) + By=o {G"^ -h) if o - 6. 

Hence, we have lifted at the multivariate level the recursive definition of 
Proposition 6 of [ABS04b]. 

Remark : The polynomial q is not positive : for G = a — 6 we have 
q{G) = - 2u' + 2v'. 
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4.2 Mixing loops and nonloops 

We now consider the polynomial Bx=y{G) := a={B(G)) where cr= is the sub- 
stitution [ya := Xa for all a &V]. A direct (static) definition is : 

B{G) = EAnB=0a;AuBw'''=(°^^[^^^l^t^"(°^^[^^^l^ 

This polynomial has two specializations : first the polynomial Q of [AvH] 
defined by Q{G, v') = t' {Bx=y{G)) where r' is the substitution : 

[u := 1; V •.= v' — 2; Xa ■= Va '■= 1 for all a € F] 

so that : 

QiG,v') = ZAHB^^iV - 2)"(GVB[AUB])_ 

Another one is the independence polynomial (Levit and Mandrescu [LM]), 
expressible by : 

IiG,v)=TjiB,=yiG)) 
where t] is the substitution [u := 0; Xa := 1 for all a GV]. 

Proposition 19 : (1) Bx=y{GVT) = Bx=y{G) for every graph G and set 
of vertices T. 

(2) A graph G without loops and the polynomial B{G) can be uniquely 
determined from p{Bx=y{G)), where p replaces w by 1. 

Proof : (1) follows from Lemma 4. 

(2) Consider two distinct vertices a and b. By looking at the ranks of the 
graphs obtained by adding loops to G[{a,b}], we see that if a — 6, then we 
have the monomials XaX^u and SxaXbU^ in p{Bx=y{G)). Otherwise, we have the 
monomials XaX^, 2xaXhU and XaX^u^-'Q 

As for Proposition 5, we ask by which algebraic transformation, one can 
recover B{G) from p{Bx=y{G)). 

Corollary 13 yields the following recursive definition : 

Proposition 20 : For every graph G : 

(1) B^ = y{0) = 1 

(2) B,^y{G) = {l + xa{u + v))B{G - a) if iV(G, a) = 0, 

(3) B^^y{G) = x,Xbii2{B,=y(G'^^ -a-b) + B^^yUG'^)'' - a - b)} 

+Xbu{B.=v{G'' -h)- B,^y{G'' -a-b)} 

+B^=y{G - a) + B^^y{G - 6) - B^^y{G -a-b) if 6 e N{G, a). 
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Wc wish to compare this definition with that presented in [AvH] for Q (and 
for graphs without loops) which we recall : 
(Ql) Q(0) = 1 

(Q2) Q{G) = u' ■ Q{G - a) if N{G, a) = 0, 

(Q3) Q{G) = Q{G -b)+ Q{G *b-b) + Q{G'''' -a) if a G N{G, b). 

Proposition 20 yields clearly (Ql) and (Q2) but, instead of (Q3) : 
(Q3') Q{G) = QiC'' -a-b)+ Q{{G''f - a - b) 
+Q{G^ - 6) - Q{G>' -a-b) 

+Q{G -a) + Q{G - 6) - Q{G -a-b) if 6 e N{G, a). 

However, Proposition 15 yields for G with a-b : 

QiG'^'' -a) = 

Q{G-a)-Q{G-a-b) + Q{G'''' -a-b)-Q{G^-a-b) + Q{{G''f-a-b), 

so that (Q3') reduces to 

Q{G) = Q{G -b) + Q{G^ -b) + QiG'''' - a). 

It remains to check that : Q{G *b-b) = Q{G^ - b). We recall that : 

G*b= {GWN{G, b))'' = G''VN{G, b). Hence : 
Q{G *b-b) = Q{G^SIN{G, b) - b) 
= Q{{G^-b)VN{G,b)) 
= Q{G'> - b) 

by Proposition 19, as was to be proved. Hence, we have established the 

recursive definition of [AvH], but not at the multivariate level. In order to 
obtain it from that of Proposition 20, we had to take u = 1 and Xa = I for all 
a. 

The advantage of the definition using (Ql), (Q2), (Q3) is that it only deals 
with loop-free graphs, whereas the definition of Proposition 20, even if used to 
compute Bx=y{G) for G without loops uses the graphs with loops (G")^ and 
G^. It proves also that Q is positive, which is not obvious from the definition. 



4.3 The independence polynomial. 



The independence polynomial is defined by 

where is the number of stable sets of cardinality k. (A looped vertex may 
belong to a stable set). Hence, we have : 

I{G,v) = n{B^=y{G)) 
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where r] is the substitution [m := 0; := 1 for all a &V]. 

We let Bi{G) = r]'{B{G)) where r]' is the substitution that replaces u by 0. 
It is a multivariate version of the independence polynomial, that can be defined 
directly by : 

whore * is the set of conditions A C V-Loops{G), B C Loops{G), GVB[A\J 
B] has no edge, so that n{GVB[A {J B]) =\ A{J B \ . From Corollary 13, we 
obtain the recursive definition 

(11) Bi{0) = 1 

(12) Bi{G) = (1 + Xav)Bi{G - a) if 7V(G, a) = 0, a is not a loop, 

(13) Bi{G) = (1 + Vav)Bi{G - a) if N{G, a) = 0, a is a loop, 

(14) Bi{G) = Bi{G -a) + Bi{G -b)- Bi{G -a-b) if 6 G N{G, a). 

However we can derive alternative reduction formulas : 

Proposition 21 : For every graph G : 
(11) Bi{0) = 1 

(15) Bi{G) = Bi{G -a) + XaV ■ Bi{G -a- N{G, a)), if a is not a loop, 

(16) Bi{G) = Bi{G -a)+ ■ Bi{G -a- N{G, a)), if a is a loop, 

(17) BiiG) = BiiG - e) - x^xtv^ ■ Bi{G - N{G, a) - iV(G, b)), 

if e is the edge a — b; (we do not delete a and 6 in G — e). 

Proof : We omit the routine verifications, which use formulas (II), (12), 
(13) and induction on size of graphs. □ 

Formulas (15), 1(6), (17) are multivariate versions of formulas given in Propo- 
sition 2.1 of the survey [LM]. 



5 Computation of interlace and other monadic 
second-order polynomials 

We consider whether and how one can evaluate for particular values of in- 
dctcrminatcs or compute (symbolically) in polynomial time the above defined 
multivariate interlace polynomials. The results of this section reformulate and 
precise results from [CMR, Mak04, Mak05, Mak06b]. In particular we consider 
multivariate polynomials with unbounded numbers of indeterminates. 

We define the size \ P \ of a polynomial P as the number of its monomi- 
als. Since monomials cannot be written in fixed size, this notion of size is a 
lower bound, not an accurate measure of the size in an actual implementation. 
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However, it is clear that the multivariate polynomial By^oiG) has exponential 
size in the number of vertices of G, and so have Bx=y{G) and B{G). Hence, 
we cannot hope for computing them in polynomial time., or by edges in other 
cases) of the considered graphs), hence as | ^ | + | _B | in the case of a mono- 
mial n ■ xaVbu^v'^- For every polynomial P{G), we denote by P{G) \ d its 
d-truncation defined as the sum of its monomials of quasi-degree at most d. 

For each rf, the polynomials i3y=o(G) |" li, Bx=y{G) \ d and B{G) \ d have 
size less than n^'^, where n is the number of vertices. Hence, asking for their com- 
putations in polynomial time has meaning. Since their monomials have integer 
coefficients bounded by n^'^, at most d occurrences of G-indexed indetcrminates 
Xa, and indeterminates u, v with exponents at most n, we can use their sizes for 
discussing polynomial time computation of these truncated polynomials. 

For their specializations P{G) where all G-indexcd indetcrminates arc re- 
placed by constants or ordinary indeterminates x,y,u,..., we have -P(G) = 
P{G) \ 0. Hence, efficient algorithms for computing d-truncations yield effi- 
cient algorithms for computing the classical (non multivariate) versions of these 
polynomials, and evaluating them for arbitrary values of x, y, u, ... 

Theorem 22 : For integers k,d, for each polynomial P among B, By^^, 
Bx^y, Bj, its d-truncation can be computed in time 0{\ V ^3d.+o{i)^ graph 
G of clique- width at most k. Polynomials g(G), Q{G) can be computed in times 
respectively 0{\ V |'') and 0(| V |^) for graphs of clique-width at most k. 

We will review clique-width in the next section. This theorem gives for 

each d a fixed parameter tractable algorithm where clique-width (but not d) 
is the parameter. See any of the books [DF] and [FG] for the theory of fixed 
parameter tractability. As a corollary one obtains the result by EUis-Monagham 
and Sarmiento [ES] that the polynomial q is computable in polynomial time for 
distance-hereditary graphs, because these graphs have clique-width at most 3, 
as proved by Golumbic and Rotics [GolRot]. 

This theorem will be proved by means of an expression of the considered 
polynomials by formulas of monadic second-order logic and for all multivariate 
polynomials representable in monadic second-order logic. 



5.1 Clique- width 

Clique-width is, like tree-width a graph complexity measure. It is defined and 
studied in [Cou97, OMR, CouOum, Oum]. 

Let G = {1, /c} to be used as a set of labels. A k-graph is a graph G given 
with a total mapping from its vertices to G, denoted by labc. We call labcix) 
the label of a vertex x. Every graph is a fc-graph, with all vertices labelled by 1. 
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For expressing its properties by logical formulas wc will handle a fc-graph 
as a tuple {V, A,pi, ...,pk) where the adjacency matrix A is treated as a binary 
relation {A{x,y) is true if A{x,y) = 1, and false if A{x,y) = 0) and p\,...,pk 
are unary relations such that Pj{x) is true iff labaix) = j. 

The operations on fc-graphs are the following ones : 

(i) For each i G C, we define constants i and for denoting isolated vertices 
labelled by i, the second one with a loop. 

(ii) For i,jGC with i ^ j, we define a unary function addij such that : 
addij{V, A, lab) = {V, A', lab) where A'{x, y) = 1 if lab{x) = i and lab{y) = j 

or vice-versa (wc want A' to be symmetric), and A'{x,y) = A{x,y) otherwise. 

This operation adds undirected edges between any two vertices, one labelled 
by i, the other by j, whenever these edges are not already in place. 

(iii) Wc let also reui^j be the unary function such that 

rerii^j {V, A, lab) = {V, A, lab') where lab'{x) = j if lab(x) = i and lab'{x) = 
lab{x) otherwise. This mapping relabels by j every vertex labelled by i. 

(iv) Finally, we use the binary operation that makes the union of disjoint 
copies of its arguments. (Hence G ® G ^ G and its number of vertices is twice 
that of G.) 

A well- formed expression t over these symbols will be called a k- expression. 
Its value is a fc-graph G = val{t). The set of vertices of val{t) is (or can 
be defined as) the set of occurrences of the constants (the symbols i and i^) 
in t. However, we will also consider that an expression t designates any graph 
isomorphic to val{t). The context specifies whether we consider concrete graphs 
or graphs up to isomorphism. 

The clique-width of a graph G, denoted by cwd( G), is the minimal k such that 
G = val{t) for some fc-expression t. A graph with at least one edge has clique- 
width at least 2. The graphs Kn, Sn-i (= Ki^n-i) have clique- width 2, for n>3. 
It is clear that clique-width does not depend on loops : cwd{GVT) = cwd{G) 
for every set of vertices T. 

The problem of determining if a graph G has clique-width at most k is NP- 

complcte if k is part of the input (Fellows et al. [Fcll+]). However, for each fc, 
there is a cubic algorithm that reports that a graph has clique-width > k ov 
produces an /(fc)-expression for some fixed function /. This latter result by 
Oum [Oum] (improved in [HO]) will fit our purposes. 

An ordered fc-graph G is a fc-graph equipped with a linear order <g on V. 
On ordered fc-graphs, we will use the variant ® of ® defined as follows : 

(iv) GIS) H is the disjoint union of G and H with a linear order that extends 
those of G and H and makes the vertices of G smaller than those of H. 
The other operations are defined in the same way. 

This extension will be used as follows (in 5.3 below) : a graph G being given 
with a fc-expression t, we replace everywhere in t the operation © by ® . The 
obtained expression t defines G and a linear ordering on its vertices. 
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5.2 Monadic second-order logic 

In a few words, monadic second-order logic is first-order logic over powersets. 
Formulas are written with special (uppercase) variables denoting subsets of the 
domains of the considered relational structures, and new atomic formulas of the 
form X & X expressing the membership of x in a set X (and for easier reading, 
also X ^ X). For more details see [Cou97] and for a use of logic in a field closely 
related to that of the present article, sec [CouOum]. 

An ordered A:-graph is handled as a relational structure {y,A,<,pi,...,pk). 
For a fc-graph, we simply omit <. Set variables will thus denote sets of vertices. 
Here are some examples of graph properties expressed in MS logic. 

That G is 3-vertex colorable (with neighbour vertices of different colors) can 

be expressed as G 1= 7, read "7 is true in the structure {V,A) representing G" 
(here <,Pi, ■■■,Pk do not matter) : where 7 is the formula : 
3Xi,X2, X3 ■ [\fx{x € XiV X G X2V X G X3)A 

Vx(-.(x e Xi A X e X2) A -.(x e A x € X3) a -.(x e Xi a x e x^)) 

A\/u,v{A{u,v) A-'{u = v) -.(u € XiAv e Xi)A^{u e X2AV € X2) 
A-.(u GX3AVG X3))]. 



That G[B] (where B CV) is not connected can be expressed by the formula 
6{X), with X as free variable : 

3r • [3x • (x e X A X e y) A 3y{y G X Ay ^ Y)A 
yx,y ■ {x G X Ay e X A A{x, y) 

=^ {{x gY Ay gY)V {x Ay ^ Y)})]. 

For B a subset of V, {G,B) 1= S{X), read "5 is true in the structure repre- 
senting G with B as value of X" iff G[B] is not connected. 

For building formulas expressing computations in GF(2), we will also use the 

set predicate Even{X) expressing that the set denoted by X has even cardinal- 
ity. This extended language, called counting modulo 2 monadic second-order 
logic is denoted by C2MS. 

Lemma 23 [CouOum] : There exists a C2MS formula p(X, Y) expressing 
that, in a graph G = {V, A) we have Y C X and the row vectors of A[X, X] 

associated with Y form a basis of the vector space defined by the row vectors of 
the matrix >1[X, X]. Hence, for each set X, the sets Y satisfying p{X,Y) have 
all the same cardinality, equal to rk{G[X]). 

Proof : We first build a basic formula X{Z,X) expressing that Z C X 
and the row vectors of A[X, X] associated with Z are linearly dependent over 
GF(2). 
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Condition Z C X is expressed by Vy ■ (y E Z => y G X). (We will then 
use C in formulas, although this relation symbol does not belong to the basic 
syntax). 

The second condition is equivalent to the fact that for each u <= X, the 
number of vertices z E Z such that A{z, u) = 1 is even. This fact is written : 

\/u-{ueX^3W- [Even{W) AWz ■ {z € W ^ z € Z A A{z, u))]). 

With X{Z, X) one expresses that Y (such that Y C X) forms a basis by : 

^A(y, X) AyZ ■ {{Y C Z A Z C X A ^{Z CY)} ^ X{Z, X)). 

We get thus the formula p{X, Y). □ 

We will say that the rank function is definable by a C2MS-formula. 

All basic results (see Section 5.4) hold for monadic second-order formulas 
written with set predicates Cardp{X) expressing that the cardinality of X is a 
multiple of p (Even is thus Card2). This extension of C2MS logic called count- 
ing (modulo) monadic second-order, can be useful for formalizing computations 
in fields GF(]3), along the lines of Lemma 23. For shortness sake we will call 
MS forrmda a formula of this extended language, (we have no reason to distin- 
guish " pure" monadic second-order formulas from those using counting modulo 
predicates), and MS logic the corresponding language. 

A function / associating a nonnegative integer f{A, B, C) with every triple 
of sets {A, B, C) is defined by an MS formula ip{X, Y, Z, U) if, for every {A, B, C) 
the number f{A, B, C) is the common cardinality of all sets D such that (G, A, 
B, C, D) N i;{X, y, Z, U). (We distinguish the variables X, F, Z, U from the sets 
A, B, C, D they denote). The generalization to functions / with k arguments is 
clear, and the defining formula has then fc + 1 free variables. 



5.3 Multivariate polynomials defined by MS formulas and 
substitutions 

For an MS formula ip with free variables among Xi, for a graph G, we 

let : 

sat(G,v3,Xi,...,X„) = {{Ai,...,An) I Ai,...,An C V, 
(G,^i,...,^„)N<p(Xi,...,X„)}. 



This is the set of all n-tuples of sets of vertices that satisfy y> in G. We can 
write it in the form of a multivariate polynomial : 

Pip{G) = Z](G,Ai,...,A„)i=¥>(Xi,...,x„) . 
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It is clear that describes exactly sat(G', (p, Xi, X„) and nothing else. 
Its set of indeterminates is Wg where W = {a;^^\ The condition 

describing the summation will be written in a shorter way fiAi, j4„). Such 
a polynomial is called a basic MS polynomial, and n is its order. 

We now show how multivariate polynomials defined by MS formulas can 
be written as specializations of basic MS-polynomials. To avoid heavy formal 
definitions, we consider a typical example : 

PiG)=j:^^A,B,C)^AyBuf(^^^'^^ (0) 

where (p{X, Y, Z) is an MS formula and / is a function on triples of sets de- 
fined by an MS formula U)- After usual summation of similar mono- 
mials (those with same indeterminates with same exponents) the general mono- 
mial of P{G) is of the form c • xaVbu^ where c is the number of sets C such 
that f{A, B, C) = p. We first observe that P{G) = <j{P'{G)) where : 

where a replaces each by 1. We are looking for an expression of P{G) as 
n{(j{Pe{G))) = fio a{Pg{G)) where n replaces each ua by u in : 

Pe{G) = J2e(A,B,CM) XAVBZcUD 

for some formula 9{X,Y,Z,U). Taking 9{X,Y,Z,U) to be A 
i}{X, Y, Z, U) would be incorrect in cases where several sets D satisfy tpiA, B, C, D) 
for a triple {A,B,C) satisfying ip. We overcome this difficulty in the following 
way : we let V be linearly ordered in an arbitrary way ; we let V'' be the formula, 
written with a new binary relation symbol denoting the order on V such that 
tp'{X, Y, Z, U) is equivalent to : 

Y, Z, U)A\/T- y, Z, T)^^^U T"] 

where U <iex T means that U is less than or equal to T in the lexicographic 
order derived from the ordering of V. This is easily expressed by an MS formula. 
The formula ijy{X,Y, Z, U) defines the function / by selecting a unique set D 
such that /(A, B,C) —\ D \. This set is unique for each linear order on V by 
the hypothesis on ip, but its cardinality does not depend on the chosen order. 
Hence we have the desired expression of P : 

P(G) = M o o'(E^(A,s,c)AV''(A,B,c,r>) xaVbZcUd)- 

These remarks motivate the following definition : 
Definition 23 a : MS-polynomials 

A multivariate MS-polynomial is a polynomial of the form : 
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where (p{Xi, ...,Xm) is an MS-formula, m' < m and f\,...,fp are MS de- 
finable functions. A multivariate MS-polynomial in normal form is defined as 
P = a o P^p where : 

^^(G):=E^(A„...,A^)<^-i^ (2) 

for ip an MS-formula and ct a substitution that can replace a variable Xa by 

1 or by an ordinary variable say u. From the above observations, it is clear 
that every MS-polynomial can be written in normal form, for graphs arbitrarily 
ordered. (The expression of P is said to be order-invariant.) 

A generalized multivariate MS-polynomial is a polynomial defined as P = 
cr o where P^ is as in (2) for ip an MS-formula and a is a specializing 
substitution. A generalized MS-polynomial may have negative coefiicients, and 
in this case, cannot have a normal form. 

Question 23 b : Is it true that every positive generalized MS-polynomial 
(i.e., that is positive for every graph) has an expression in normal form ? 

We say that P is of order m, if it can be expressed as P = cr o where ip has 
TO free variables. Hence a polynomial of the form (1) above is of order m + p. 

Then for a graph G with n vertices and P defined by (1), P(G) has size 
at most 2™ degree at most n{m' -\- p), positive coefficicints of value at most 
2™". These bounds will be useful for evaluating the cost of computations of 
truncations of polynomials. 

Transformations of MS-polynomials. 

We now show how some specializations can be reflected by transformations 
of the defining formulas. We review some cases which arised in the present 
article, by taking a polynomial P of the form (0). For each case 1,...,4 we 
denote by Pi the polynomial obtained from P. 

Case 1 : Xa 0. We have : 

where ip'{Y, Z) is defined as 3X ■ [ip{X, Y, Z) hMz ■ z ^ X\ so that >p'{B, C) is 
equivalent to (^(0, B, C). The function /(0, Y, Z) is defined by 3X-[^{X, Y, Z, U)A 
Vz-z^ X]. 

Case 2 : u := 0. We have 

^2(G) = Ev'(A,B,C) ^AVB 

where ip' {X, Y, Z) is </p(X, F, Z) h3U ■ [V'(X, Y,Z,U) A^z ■ z ^U], which is 
equivalent to p{X, Y, Z) A " f{X, Y, Z) = 0" . 
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Case 3 : Xa := 1. We have 

The sets A are, hke the sets C, "invisible" in the monomials. However, 
they play a role. They contribute to the multiplicity of the monomials ysu^, 
for p = f{A,B,C). (This case has been considered above at the beginning of 
Section 5.3). 

Case 4 '■ Xa ■— Ua- Here 

We assume here, and this is enough for this article, that Lp{A, B, C) implies 
that A and B are disjoint. Then to reach the general syntax we write P4{G) as 
follows 

where D stands for Au B, ip'{W,Y,Z) is chosen to be equivalent to the 
formula Y CW A (p{W — Y, Y, Z) and the function g is defined by : 

g{D,B,C) = f{D-B,B,C) 

whence also by a formula i^'{W, Y, Z, U) equivalent to : 
3X • [il){X, Y,Z,U) AX CW A''Y = W - X"]. 

Lemma 24 : Each polynomial P among B, By^o, B,.=y, Bj is an MS- 
polynomial. For every graph G, ordered in an arbitrary way, we have an expres- 
sion of P in normal form P{G) = a o Pg(G) for an MS formula 6 expressing 
properties of ordered graphs, and some specializing substitution a. 

Note that P{G) = cF{Pe{G)) for every linear order on G. 

Proof : We only consider B. The other cases follow by the techniques 
presented above. We recall the definition of B: 

We let 

PoiG) = 'Ee(A,B,C,D) xaVbUcVd 

where e{A, B, C, D) holds Af^ B = 0,G C A\J B,D = A\J B - C, G \s 
the smallest basis of the vector space spanned by Mgvb[A U B, ^ U S], where 
Mgvb is the adjacency matrix of GVB. By "smallest" we mean with respect 
to the lexicographic ordering derived from the ordering of G. It follows that 
I G 1= rk{GVB[A U B]) and | D \= n{GVB[A U B]). By Lemma 23, one can 
express these conditions by an MS formula 0. 

Hence, B = a o Pg where a replace Ua by u and Va by 
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5.4 The Fefermann and Vaught paradigm applied to MS- 
polynomials 

We now show how multivariate polynomials defined by MS formulas can be com- 
puted "efficiently", not by reduction formulas, because they yield (when they 
exist) exponential computations in general, but by induction on ^-expressions 
defining the considered graphs. This will only apply to graphs of clique-width 
bounded by a fixed value. We will use the Fefermann- Vaught paradigm, pre- 
sented in detail by Makowsky [Mak04, Mak05]. 

We need operations that manipulate sets of q-tuples, in particular, those of 
the form sat(G, tp, Xi, Xq), and, equivalently as we will see, the polynomials 
PAG)- 

For sets R, S and S' C P{V)'', we write 

R = S W S' a R = S U S' and S n S' = 0, and 

R = S^S' use V{Vi)i, S' C P(F2)«, Vir\V2 = 0, and R is the set of 
g-tuples (^1 U Bi, Ag U Bg) such that {Ai, Ag) e S and (Bi, Bg) e S'. 

For each S C 'P(y)', we let P{S) be the multivariate polynomial : 

P{S) = Y^(Ai,...,Aq)&S ^aI---^aI ■ 

Clearly, P^(G) = P(sat(G, >fi,Xi, ...,Xq)). The following is clear : 

Fact 25 : For 5, S" C V{V)i, we have P{S tb) S') = P{S) + P{S') and 
P{SmS')=P{S)-P{S'). 

We denote by Ufe the (finite) set of imary operations allowed in fc-expressions. 
We denote by MS{k,q) the set of MS formulas written with the basic symbols 
=,G,^,Cardp and the relation symbols A, <,pi, ...,pk (hence able to express 
properties of ordered fc-graphs) with free variables in the set {Xi, Xq}. The 
following theorem discussed in [Mak04] is proved in [Cou97, CouMos] in closely 
related forms : 

Theorem 26 : For every k, q, for every formula ^ in MS{k, q), there exists a 
finite subset $ of MS{k, q) containing ^ and satisfying the following properties: 

(1) For every ip G ^ for every op G Ufe there exists a formula (p°P G $ such 
that, for every ordered fc-graph G : 

sat(op(G), if, Xu Xg) = sat(G, ^"p, Xu Xq). 

(2) For every ip G ^ there exist p and ...,6p,ipi, ...,tjjp) G such that 
for disjoint ordered fc-graphs G and H : 

sat{G®H,<p,Xi,...,Xq) = 

Wi<i<pSat(G, ei,Xi, Xg) m sat{H, ij^^Xi, Xg). 

These statements also hold for (unordered) fc-graphs and the operation ® 
instead of © . 
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Let $ be a set of formulas as in Theorem 26. We get a finite family of poly- 
nomials (P,p),^g$ that satisfy mutually recursive computations rules. Actually, 
the recursive rules apply to the family of polynomials {a o P^)^^^ where cr is a 
specializing substitution. We recall that by a polynomial wc mean (ambigously) 
a mapping P associating with a graph G a polynomial in Z[[/ U Wg]- 

Corollary 27 : Let $ satisfy the properties of Theorem 26. Let cr be a 
specializing substitution. We have the following computation rules : 

(1) For every y G for every op € Ufc for every ordered fc-graph G : 

[a o P^){op{G)) = (ctoP^op)(G). 

(2) For every <^ G $ for every disjoint ordered fc-graphs G and H : 

(<7 o P^)iG®H) = Ei<i<p(^ ° PoMG) ■ (a o P^,){H). 
where (p"P and {6i, ...,9p,ipi, ipp) are as in Theorem 26. 

Proof : If (7 is the identity substitution, then (1) and (2) are direct trans- 
lations of (1) and (2) of Theorem 26. 

Since ao{P+Q) = aoP+aoQ i.e., {aoP+aoQ){G) = a{P{G))+a{Q{G)) = 
a{{P+Q){G)) and similarly, ao[P-Q) = (croP) • (cjoQ), for every substitution 
a and every two polynomials P, Q, the equalities extend to the general case as 
stated. □ 

Hence, this corollary concerns all multivariate polynomials described in 5.3. 
We will use it for their computation. 



5.5 Computing polynomials in polynomial time 

We discuss the computation of polynomials written in the form a o P^ for P^ 
a basic MS-polynomial of order q, and a substitution a. This will also apply 
to evaluations of polynomials where all indeterminates are given some numeric 
value, either integer, real or complex. 

The number n will denote the number of vertices of the considered graph 
G. In evaluating the cost of a computation as 0(e) where e is a nonnegative 
expression, we omit the case where e might have value (rigourously, we should 
write 0{e + 1)). Similarly log(e) stands for Maa;{l, log(| e | +1)}. 

Let a graph G be given by a fc-cxprcssion t. Wc "order" t into t which 
defines G with a linear order of its vortices (cf. Section 5.1). 

It is clear that for each constant, i or i^, each polynomial (a o P^)(i^) can 
be computed from the definitions. We can thus compute for each subterm s of 
t the family of polynomials {{a o P^p){val{s)))^p(z^. In particular, at the end 
of the computation, one gets (tr o P^){G) for all e (Which is too much. A 
method for restricting similar computations to their necessary parts is described 
in [CouMos] or in [Mak04], Definition 4.17.) 

This computation uses at most n- | $ | times the computation rules of 
Corollary 27(2), because in a term, the number of occurrence of ® is s — 1, 
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where s is the number of occurrenees of constants, which is equal to | y |= n. 
Here, | $ | is constant. The computation time is bounded by 2n-CG-Pmax- I ^ I) 
where Pmax is the maximum value of p in the rules of Theorem 26(2), and cg 
bounds the cost of the addition and of the muhiphcation of two polynomials. 
This bound depends on G. When do we obtain a polynomial algorithm in n? 

We first precise the way we will count the cost of operations on polynomials. 
Using unit cost measure, we will count for one the cost of each basic operation 
on numbers : comparison, addition, substraction, multiplication. The cost of 
evaluating .t™ for a positive integer m is thus 0(log(m)). In the computation 
of a truncated polynomial P \ d, we will consider d as "relatively small" and 
fixed, like | $ | (and actually much smaller in potential applications.) Hence 
we will count for one the cost of computing .t™ for m < d. 

However, for computing a polynomial or evaluating it for integer values of 
the arguments, we can also use the real cost measure and consider that the 
cost of a comparison, an addition or a substraction of two positive integers x 
and y is 0(log(| x | + | y |)) and that of their multiphcation is 0{\og{xy)). 
Coefficients of polynomials may be exponentials in the sizes of the considered 
graphs. However, in situations where the absolute values of coefficients and 
exponents are no larger than 2p'"-' for fixed polynomials p, a polynomial bound 
on computation time with respect to unit cost measure remains polynomial with 
respect to real cost measure. The exponents of the polynomial bounds are just 
larger. 

We will base the following estimations of the cost of computations on straight- 
forward data structures : a polynomial is a uniquely defined list of monomials 
sorted by increasing order of quasi-degree, where two monomials of same quasi- 
degree are ordered lexicographically. Each monomial is written in a canonical 
way by means of a fixed ordering of indeterminates. We deal with monomi- 
als with a variable number of indeterminates, however, this number is always 
bounded by n- \ X \ + \ U \ and where is the set of G-indexed indetermi- 
nates, and U the set of others. 

The basic operations on pairs of monomials m, m' are comparison, summa- 
tion of coefficients if m, m' are similar monomials (if they have same indetermi- 
nates with same respective degrees), and multiplication. For a monomial m, we 
denote by v(m) the number of its indeterminates. The costs of these operations 
are respectively 0(v(to) -1- v(m')), 1, 0(v(m) + v(m')).We denote by v(P) the 
number of indeterminates in a polynomial P. 

Lemma 28 : (1) For every P, Q, d, if v(P), v{Q) < v^ax we have: 
{P + Q)\d^Pld + Q\'d and (PQ) \ d = {{P \ d) ■ (Q \ d)) \ d. 

(2) Computing P + Q takes time 0(vmax • (I P I + 1 Q I))- 

(3) Computing PQ takes time : 0(vmax- | P ^ • | Q |) if | P |<| Q | . 

Proof : (1) Clear from the definitions. 
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(2) Note that |P + (3|<|P| + |Q|. The addition of P and Q is done 
by interleaving their hsts of monomials and by adding the coefficients of similar 
monomials. This gives the result by the remaks on the cost of operations on 

monomials. 

(3) Let I P |<| I and v(P),v(Q) < v^^^. Note that | PQ |<| P | • | Q |. 
We compute PQ by multiplying | P | times the polynomial Q by a monomial 
of P, and by performing | P | — 1 additions of polynomials of size at most 

I P I • I Q |. The time taken is at most 

0{W^^- I P I • I I +(| P I -1) • V^a.- I P I • I I) = 0{V^^- I P |2 . I Q I). □ 

Remark : If in (3) P and Q are positive, one gets the bound 

0(vmax- I P I • I PQ I) because all intermediate results have size at most 

\PQ\- 

Theorem 29 : Let k, d be fixed integers. The rf-truncation of a generalized 
multivariate MS-polynomial P{G) can be computed in time 0{n^'^~^'^^^^) for 
every graph G of clique-width at most k. 

The constants hidden by the O-notation depend on P and k. The particular 
case of evaluations (for numerical values of indeterminates) will be discussed 
later. Closely related formulations of this theorem are in [CMR, Mak04, Mak05] . 
For a polynomial P, we define |1P|| := (| P |, dcg(P), Cmax(P)) where dcg(P) 
is the degree of P and Ci„ax(P) the maximum absolute value of its coefficients. 
Triples of integers are ordered componentwise. 

Lemma 30 : Let P be a polynomial and cr be a substitution such that 
||f(a;)|| < {smax,dmax,Cmax) for every indeterminate x. The polynomial a o P 
satisfies : 

||aoP|| < (I P I •(s„„,)'^«s(^),deg(P)-d„„,,| P I •C^a.(P)-(s™axW)d<'s(P)). 
Proof : Easy verification. □ 

It follows from previous observations (sec the remark after Question 23 b) 
that if Pip has order m, then for every graph G, the polynomial a{P^{G)) has 
size and coefficients bounded by 2*^'"^ and degree bounded by n{ni'+p)-dmax = 
0{n) [m' < m and p are as in that remark). We can thus use the unit cost 
measure. 

Proof of Theorem 29 : Let fc, d be integers, let P be a polynomial ex- 
pressed as cr o for an MS formula ^ and a substitution cr. We aim at 
computing its d-truncation. Let $ be the corresponding set of formulas as in 
Theorem 26. We observe that by Lemma 28 (1), Corollary 27 yields : 
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(1) For every £ $, for every op € XJk for every ordered fc-graph G : 

a{P^{op{G))) \ d = a{P^a,{G)) \ d. 

(2) For every e $ for every disjoint ordered fc-graphs G and H : 

a{P^{G®H)) \ d = Ei<i<piHPeAG)) \ d) ■ {a{P^AH)) \ d)) \ d. 

Note that we do not have : a{P{G)) \ d — a{P{G) \ d) in general. As 
observed above, the time to compute a{P{G)) is 2n ■ cq • Pmax- I I, where 
cg bounds the costs of adding and muhiplying the polynomials occuring in 
recursion rules (1) and (2). We need only count multiplications which are more 
costly than additions and in proportional number. Let us assume that aU the 
polynomials a{Pip{G)) arc in Z[Xg U U]. By Lemma 28 we have, for the d- 
truncations of all such polynomials P, Q : 

(i) CG = 0(v„,ax- \P \d\^ ■\Q\d\) with 

(ii) v„,ax = n-\X\ + \U \= 0(n), and 

(iii) I P rd |=0(n2'^+lt^l). 

For proving (iii), we note that a monomial of P f d is a product of at most d 
factors of the form a;^ and of at most | U \ factors of the form , in both cases for 
s < dcg(P) -dmaa;. There are n- | X | • deg(P) -dmax factors of the form .t^ and 
deg(P)-dmax factors for each m. Hence \ P \d\= 0((n-deg(P))'^-deg(P)l^l) 
= 0(n2'^+l^l) since deg(P) = 0(n) . 

This gives for 2n • cq ■ Pmax- | $ | the bound 0{n'^+^'^+^\^\). 

We must take into account the cost of building for a graph G of clique- 
width at most k a clique-width expression. In cubic time, one can construct 
for graphs of clique- width at most k an /(A;)-expression, for a fixed function 
/ ([HO, Oum]). This suffices for our purposes. The total time is thus 0(n*) 
where t = Maa;{3, 2 + 6d + 3- | C/ |}. (This bound applies iid=\U \= 0.) □ 

We extend this result to numerical evaluations. Let P{G) be a polynomial in 
Z[Xg U U] to be evaluated for a graph G. An evaluating substitution v replaces 
indctcrminates by integer, real or complex values. Let v be such a mapping. 
One can consider v(P{G)) as a polynomial reduced to a constant, that is the 
desired value of P{G) for the values of indeterminates specified by v. Note 
that iy(P{H)) is well-defined for every graph H with set of vertices included in 
the set V of vertices of G. This remark will be useful for the computation of 
u{P{G)) by induction on the structure of G using Corollary 27. The costs of 
computations are the same for polynomials with integer, real or complex values 
since we use unit cost measure. 

Corollary 31 : Let k be an integer. For every generalized multivariate 
MS-polynomial P and every evaluating substitution, the corresponding value of 

P{G) for a graph G of clique-width at most k can be computed in cubic time 
in the number of vertices of G. It can be computed in linear time if the graph 
is given by a fc-cxprcssion. 

Proof : Let v be an evaluating substitution. It associates a number with 
each u'mU and each Xa in Xg, for a given graph G. As in the proof of Theorem 
29 we have : 



31 



(1) For every e for every op S Ufe for every ordered fc-graph H with 

u{P^{op{H))) = u{P^o,{H)). 

(2) For every G $ and all disjoint ordered fc-graphs H and H' with 

Vh,Vh' C Vg : 

Here wc compute vahics, not polynomials. The cost is thus 2n -Pmax' | ^ | = 
0{n) assuming known a fc-expression defining G. □ 

We now precise the bounds for the polynomial B and its specializations. 
Theorem 22 is actually a corollary of Theorem 29. 
Proof of Theorem 22 : 

For every graph G with n vertices, we have the following bounds : 

\\B{G)\\ < (3", 2n, 1), I B(G) \ d \= 0(n'^+2), 
||B^=o(G)|| < (2», 2n, 1). I S,=o(G) \ d \= 0(n^+i), 
||B,=j,(G)|| < (n2",2n,2"),| S,=y(G) f d |= 0(n'^+2), 
\\Bi{G)\\ < (2", 2n,l), I i?7(G) f d |= ©(n-^+i), 
lig(G)|| < (n2,2n,2"), 
\\Q{G)\\ < (n+l,n,3"). 

As in the proof of Theorem 29, we need only bound the costs cg of the multi- 
plications of polynomials. From these evaluations, wc gc^t the bounds 0{n^'^^^) 
for B and B^^y,0{n^'^+^), for By^o and Bj, O(n^) for q and Oin^) for Q 
for the last two, Vmax is constant. □ 

Remark about the size of constants. 

Sets $ in Theorem 26 are very large if they are constructed in a blind 
manner from an MS formula : the size is a tower of exponentials proportional 
to the quantification depth of formula ^. However, if alternatively a family of 
polynomials a o satisfying Corollary 27 is constructed directly, by using our 
knowledge of the meaning of the properties defined by formula ^, then one may 
obtain a usable recursive definition. Hence the above estimations leave a great 
space for improvements. 



6 Conclusion 

We have defined a multivariate interlace polynomial that generalizes the existing 
interlace polynomials. The multivariate methodology puts in light the mean- 
ing of polynomials. Classical polynomials are degraded versions of multivariate 
ones. The multivariate approach is well-adapted to the logical description of 
polynomials. And the use of monadic second-order logic yields FPT algorithms 
for evaluating polynomials at particular values or for computing significant por- 
tions of them, called truncations. 
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For computing such polynomials in full, one might use linear delay enumer- 
ation algorithms and try to obtain monomials one by one, by increasing degree, 
with a delay between two outputs linear in the size of the next output. Such 
algorithms arc considered in [Cou06] and [Bag] for MS definable problems and 
graphs of bounded clique- width. 

Another research perspective consists in enriching the notion of configura- 
tion. In this article, in particular in Section 5, a configuration is an m-tuple of 
subsets for fixed m. One could try to extend the methodology of Section 5.3 to 
more complex configurations like partitions of unbounded size or permutations. 

Finally, in order to build a zoology as opposed to maintaining a zoo (as 
written by J. Makowsky [Mak06a]) it is important to relate the various poly- 
nomials by means of algebraic reductions (specializations), or logical reductions 
or transformations of other kind. 
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8 Appendix : Tutte polynomial 

In this section we apply our logical tools to the Tutte polynomial of matroids. 
We define a multivariate Tutte polynomial "better" than Sokal's. We first es- 
tablish a lemma. 
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8.1 Enumerating polynomials 

The enumerating polynomial of a set S C V{V) is defined as Enum{S) = 
^Aes ^® write Enumx{S) to specify the generic type x of indeterminates. 
(Hence Enumy{S) = Y^^agsVa)- 

For a polynomial P, wc let P~ (rcsp. P"*") denote the polynomial obtained 
by replacing each indeterminate x (of any type) by a; — 1 (resp. by a; + 1). 

If S C y, we let Sub{B) — {AlAC B}. 

Lemma A.l : For every S C VCV), the following are equivalent : 

(1) 5 = Sub{B) for some B CV, 

(2) Enum{S)~ is positive, 

(3) Enum{S)~ = xb for some B CV. 

Proof : (!)=> (3) =^(2) is straightforward. The set B in (3) is the same 
as in (1). 

(3) =^(1) : We note that Enum{S) = {Enum{S)-)+ ; ii Enum{S)- = Xb 
then Enum{S) = Xb^ = Enum{Suh{B)) hence S = Suh{B). 
(2)^> (3). By induction on the cardinality of V. 

Wc let 5 \ a denote {A - a \ a & A e S) and S - a = {A \ a ^ A & S). We 
have thus 

Enum{S) = Xa ■ Enum{S \ a) + Enum{S — a), and 
Enum{S)~ = {xa — 1) • Enum{S \ a)~ + Enum.{S — a)~ . 
It is positive by hypothesis. So is Enum{S \a)~ (just look at monomials 
with Xa), hence S\a = Sub{B') for some B' CV — a. 

Then Enum{S — a)~ — Enum{S \a)~ is also positive. Since this is the case 
for Enum{S \ a)~ the same holds for Enum{S — a)~ . 
Hence S — a = Sub{C) for some C CV — a. Hence 

Enum{S)~ = {xa — +xc- That Enum,{S)~ is positive implies B' = C. 
Hence S = Sub{B) where B = B' U {a} and Enum{S)- = ccs .□ 



8.2 The Tutte polynomial of a matroid 

We represent the Tutto polynomial of a matroid as a specialization of another 
multivariate MS-polynomial than the one defined by Sokal in [Sok] recalled in 
the introduction. We establish at the multivariate level that its coeflScients are 
positive, with help of the notion of an active element with respect to a basis 
and a linear ordering. 

We consider the Tutte polynomial for matroids but the application to graphs 
is immediate. Our reference is the book edited by N. White [Whi]. Wc refer to 
this book for the basic definitions on matroids and some results. A matroid M = 
{E, I) with base set E (we can think of E as the edge set of a graph) is considered 
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as a relational structure where / is a set predicate defining the independent sets. 
MS logic can be used to describe matroids (see [Hli03], [CouOum]). For instance, 
the bases of M are characterized by the MS formula P{X) : 

I{X) A VZ • [{X C Z A I{Z)} ^ Z C X]. 

The rank polynomial of M is defined as : 

R{M) = 

where r(M) is the rank of M, r{A) is that of A and n{A) :=| A \ -r{A). 
The Tutte polynomial is defined as : 

T(M) := R{M)- := Y.a^e{x ' If^^^-'^^^Xy - l)"^^) • 

Whereas R{M) is clearly positive, this is not obvious for T{M) from this 
definition. 

We let (T be the substitution that replaces each Xa by x and each ya by y. 
In our logical setting, the rank polynomial can be expressed as <j o i?, for M 
linearly ordered in an arbitrary way where : 

RiM) = 'E^(A,C\D)^CyD, 

and the MS-formula <f{A, C, D) expresses the following conditions : 

(a) A<ZE, 

(b) there exists a set Z such that Z C 7I, Z is a maximal independent subset 
of A, that is lexicographically minimal with these properties, 

(c) C C — A is lexicographically minimal such that Z U C is a base of M, 

and 

(d) D = A-Z. 

Hence (C, D, Z) is associated in a unique way with A and | C |= r{M)—r{A), 
\D\=\A \ -r{A) = n{A). It follows that R{M) = a{R{M)). 

Hence, the rank polynomial is an MS-polynomial. For proving that the 

Tiitte polynomial is positive, we will replace R by another multivariate poly- 
nomial, in such a way that we can deduce from Lemma A.l that the Tutte 
polynomial is positive. Our main tool is the following proposition from [Bjo]. 
We need some definitions relative to an ordered matroid M = {E, I, <) (i.e., E 
is linearly ordered). 

For every base B, an element e of E — B is externally active with respect to B 
if it is the least element with respect to < of its fundamental cycle (the unique 
cycle included in i? U {e}). Their set is denoted by EA{B). Dually, an element 
6 of -B is internally active in B if it is the least element of its fundamental 
cocycle. Their set is denoted by IA{B). 

Proposition A. 2 ([Bjo], Proposition 7.3.6) : For every ordered matroid 
M={E,I,<) : 
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V{E) = ^f3(B){A I B - IA{B) CACBU EA{B)}. 



Wc recall that W denotes a disjoint union and that (3{B) expresses that B 
is a basis. Hence every subset A of E can be decomposed in a unique way as : 

(i) A={AnB) U{A-B) 

(ii) with ADB D B - IA{B) and A- B <Z EA{B), for a unique base B. 

(iii) We IctC = B-A. Hence C C IA{B). Then r{A) ^\AnB\ .(See [Bjo] 
for the proof). It follows that | C |= r(M) - r{A). 

(iv) Let D = A - B: then n{A) =\ D \. 

The niouomial xcUd satisfies : 

(v) =a;'^(^)-'^(^)j/"(^). 
We define : 

^{M) = E/3(B) T,e{B,c,D) ^cyo, 

where 6'(S, C, L>) expresses that C C IA{B) and £> C EA{B). 
Claim : i?(M) = a{R{M)). 

Proof : Let A C £ and m{M,A) = 2;'-(M)-r(A)yn(A) ^j^g corresponding 
monomial in the expression of R{M) as X^acb ^{M, A). 

Let _B be the unique basis satisfying conditions (i)-(ii). Let C and D be 
associated with A by (iii) and (iv). Then, by (v), m{M, A) is the image by a of 
a unique monomial of R{M). 

Conversely, let xcUd be a monomial of R{M). Let A = {B — C) L) D. Con- 
ditions (i)-(iv) hold, hence {B,C,D) is the unique triple associated with A by 
these conditions, and thus crixcyo) = ni{M,A). This proves the claim. □ 

Note that for each basis _B,the condition 0{B, C, D) expresses that C and 
D can be chosen independently, with C C IA{B) and D C EA{B). Hence 
R{M) = a{R{M)) where : 

R{M) = Y.f){B) Enum^{Sub{IA{B))) ■ Enumy{Sub{EA{B)). 
Let us define 

T{M) = E/3(S) XiA(B)yEA(B)- 

This is our alternative multivariate Tutte polynomial for an ordered matroid 
M. Clearly, 

R{M)- = 

= {E/3(s) Enum^{Suh{IA{B))) ■ Enumy{Suh{EA{B))}- 

= E/3(B) Enum^{Suh{IA{B)))- ■ Enumy{Sub{EA{B))- 

= T,f3{B) Enumx{{IA{B)}) ■ Enumy{{EA{B)}) = E/3(b) xiA(B)yEA{B) 
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= T{M). 

We have : _ _ _ 

RiM)- = a{R{M))- = a{R{M)-) = a{T{M)) = T{M). 

Hence T{M) and T{M) are both positive. We have Ufted at the multivariate 
level, the well-known identity : 

= R{M)- 
because we have T{M) = R{M)- . 

Computing and evaluating the Tutte polynomial. 

The polynomials 

R{M) = T,/3iB) Ee(B,c,r>) ^cVd, and f (M) = T,/3{b) XiA(B)yEA(B) 

are thus MS-polynomials, because the notion of a basis, the sets IA{B) 
and EA{B), the conditions 6{B, C, D) are MS-expressible. So are R{M) and 
T{M). Some consequences regarding their computations follow along the lines 
of Theorem 29. 

This approach is developped by Makowsky in [Mak05] for graphs of bounded 
tree-width. (The application of Theorem 29 is not immediate because a linear 
ordering of the edge set must be incorporated to the relational structure in such 
a way it remains of bounded tree-width.) Polynomial algorithms that compute 
the Tutte polynomial of graphs of bounded tree- width are given in [And, Nob]; 
these algorithms do not use monadic second-order logic. However they use so- 
called "splitting" formulas which can be seen as particular cases of equalities 
(1) and (2) of Theorem 26, which are central for the proof of Theorem 29 (as 
pointed out in [Mak05]). 

Gimenez et al. give a non-polynomial algorithm that computes the Tutte 
polynomial of graphs of bounded clique- width ([GHN]). Theorem 29 is not ap- 
plicable in this case because edge sot quantifications are necessary and MS2 logic 
does not fit with bounded cHque-width. 

For matroids of finite branch- width represented on finite fields, Hlineny gives 
a polynomial algorithm [Hli06] that follows the ideas of those of [And, Nob]. 
However MS properties of such matroids can be checked in polynomial time 
(see [Hli03]), hence the method of Theorem 29 may perhaps be applied also to 
these matroids. 
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